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Abstract The first near-side X-class fiare of ttie Solar Cycle 24 occurred in 
February 2011 and produced a very strong seismic response in the photosphere. 
One sunquake was reported bv lKosovichevl f Astrophys. J. Lett. 734, L15 J2011]) , 
follo wed by the discovery of a second sunquake by Zharkov, Green, Matthews 
et al. ( Astrophys. J. Lett. 741. L35 J2Qlip . The flare had a two-ribbon structure 
and was associated with a flux rope eruption and a halo coronal mass ejection 
(CME) as reported in the CACTus catalogue. Following the discovery of the 
second sunquake and the spatial association of both sources with the locations 



of the feet of the erupting flux rope (Zharkov, Green, Matthews et al. Astrophys. 



J. Lett. 741, L35. I201T]) . we present here a more detailed analysis of the observed 
photospheric changes in and around the seismic sources. These sunquakes are 
quite unusual, taking place early in the impulsive stage of the flare, with the 
seismic sources showing little hard X-ray (HXR) emission, and strongest X-ray 
emission sources located in the flare ribbons. We present a directional time- 
distance diagram computed for the second source, which clearly shows a ridge 
corresponding to the travelling acoustic wave packet and flnd that the quake at 
the second source happened about 45 seconds to one minute earlier than the 
flrst source. Using acoustic holography we report different frequency responses 
of the two sources. We flnd strong downflows at both seismic locations and a 
supersonic horizontal motion at the second site of acoustic wave excitation. 
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1. Introduction 

Sunquakes are observed as photospheric ripples, which accelerate radially out- 
ward from a source region. The theoretical prediction that sunquakes should 



be produced by the energy released during major solar flares ( [Wolff, 1972 ) was 
supported by their discovery on the Sun bv 'Kos ovichev and Zharkoval (jl998p . 
The acoustic nature of quakes has been well established since their discovery. 
However, the exact physical mechanism behind their excitation is still debated 
with several theories currently under consideration. Observations of sunquakes 
are relatively rare, possibly due to the difficulties of detecting the photospheric 
ripples, and helioseismic methods such as time-distance diagram analysis and 
acoustic holography are employed to look for evidence of acoustic emission. 
With only a small number of such events verified so far (see Besliu-Ionescu, 
Donea, Gaily et a/. n2QQ5l [Donea, Besliu-Ionescu, Gaily et oTf I2QQ6I for some 
examples of known quakes from the last solar cycle) the new solar cycle and the 
virtually c ontinuous high-resolution data of SDO/HMI ( Scherrer, Schou, Bush 
et a/., 2Q12; |Schou, Scherrer, Bush et al.^ 2Q12D mean that sunquake detections 
should increase in the coming years. 

The 15 February 2011, X2.2 class fiare occurred in NO A A active region 11158 
and was the first in the much delayed rising activity phase of the new Solar Gycle 
24. The active region started emerging in the eastern hemisphere on 10 February 
2011, with two bipoles emerging side by side creating a complex multipolar 
region. As the active region evolved through both emergence and cancellation 
events, the coronal loops became increasingly sheared resulting in a number of 
G-class and M-class fiares occurring from 13 February onward, culminating in 
the X-class event with GOES fiux peaking around 01:55 UT on 15 February. 
The X-class fiare was a long-duration fiare with an impulsive phase, as observed 
in GOES 1.0 to 0.8 A soft X-ray data, lasting from 01:46 to 01:56 UT, and 
integrated HXR emission observed by RHESSI up to ^ 100 keV, peaking just 
before 01:55 UT. The strongest HXR emission was produced in the 6 — 25 keV 
energy band. 

In the photosphere, the fiare exhibited a classic two-ribbon pattern, which 
is most clearly seen in the SDO/HMI line-of-sight magnetic field and ve- 
locity running-difference images as well as Hinode/SOT Ga ll K observa- 
tions {e.g. Figure [T]). Though less pronounced, the ribbons are also present 
in the SDO/HMI continuum data. Spatially, hard X-ray emission was situ- 
ated primarily along the ribbons. The halo GME associated with this ffare 
was detected by GAGTus software and is listed in the LAS GO catalogue 
(http://sidc.oma.be/cactus/catalog/LASCO/2^_0/ qkl/2011/02/). 

The ffare produced a strong seismic response ( [Kosovichev, 201 1"! , with rip- 
ples travelling outward from the source clearly seen in HMI velocity difference 
data (see, for instance, the online movie in the above article). Using acoustic 
holography [Zharkov, Green, Matthews et al.\ (|201ip have shown that a second. 



apparently weaker, source of acoustic waves is present and have shown that the 
sunquakes occurred at the foot-points of a ffux rope. The presence and location 
of the ffux rope has been deduced based on a straight-forward observational case, 
making use of the photospheric magnetic ffux distribution, sigmoidal structure. 
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Figure 1. Sunquake source locations determined from egression power computed at 6 mHz 
and time-distance overplotted on egression power snapshot (a), continuum intensity (b), mag- 
net ogram (c) and velocity (d) images. Panel (e) shows the flare-induced changes in magnetic 
fleld computed as the difference between magnetic data in panel (c) and the 20-minute average 
of HMI line-of-sight magnetorgrams before the flare onset. The orange contours in (a) - (e) 
correspond to 01:51:27 6-mHz egression power contours at 2.5 and 3 times the quiet-Sun 
egression at this frequency. The red stars in all images mark the locations used for computing 
time-distance diagrams. Panel (f) is Ca ll K taken by Hinode/SOT at 01:51:39 UT with HXR 
contours deduced from RHESSI data in 25 - 50 keV range overplotted. The HMI data in panels 
(a)-(e) are Postel projected so the distances along x- and 7;-axes are plotted in megameters. 
Arcsecond coordinates are given along the axes in panel (f). 

chromospheric and coronal changes during the eruption. The observational inter- 
pretation of the presence of a flux rope is supported by non-hnear magnetic-field 
mod elling from HMI vector magnetogram data ( Schrijver, Aulanier, Title et al.^ 
2011; |Sun, Hoeksema, Liu et al, 2Q12D . 

In this article we analyse seismic measurements and report the properties of 
the detected seismic sources and consider associated changes in photosphere. The 
data and methods are described in Section [2j results are presented in Section [3l 
followed by discussion and conclusions. 
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Figure 2. Time-distance diagrams for both seismic sources. The first two plots from the 
left correspond to the first (eastern, strong) source, the following are the western source. 
Theoretical time-distance curves are overplotted in white. Locations of the sources are marked 
as red stars in Figure ^ 

2. Data and Methods 

We use full disk SDO/HMI intensity, dopplergram, and line-of-sight magnetic 
field data at 45-second cadence to produce three hour-long datacubes, which 
are extracted by remapping and de-rotating the region of interest using Postel 
projection and the Snodgrass differential rotation rate. The spatial resolution is 
0.04 degrees per pixel. The hard X-ray data come from RHESSI (Lin, Dennis, 
Hurford et a/., 2002), which observed the flare from the pre-cursor phase begin- 
ning at 01:27 UT until 02:30 UT, covering the entire impulsive phase. We used 
the CLEAN algorithm to produce images at between 20- and 40-second cadence 
covering the duration of the flare. 

We apply acoustic holography to calculate the egression power maps from 
observations . The holography method (|Braun and Lindsey[ [T999| D onea, Braun, 
and Lindsey"E9991 |Braun and Lindsey| [20001 [Lindsey and Braun| l2QQQl [2QQTj) 
works by using Green's function [G+(|r — r'|, t — t')] which prescribes the acoustic 
wave propagation from a point source, to essentially "backtrack" the observed 
surface signal [?/;(r,t)]. This allows us to reconstruct egression images showing 
the subsurface acoustic sources and sinks: 

H+{v,z,t)= Ut' j dVG+(|r-r'|,t-tW,0, (1) 

where a, h define the holographic pupil. The egression power is then defined as 

P(^,r) = j |i7+(r,^,t)pdt, (2) 

and can be viewed as a proxy for acoustic energy emitted from a location around 
{z^y) over the integration period. In the above, r represents horizontal position, 
z depth, and t time. 

In this work, the egression power is computed for each integral frequency from 
3 to 10 mHz, by applying 2-mIIz frequency bandwidth filters to the data (for 
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Figure 3. Source 1 (left) and 2 (right) velocity transients associated with the quake. Top 
row: HMI velocity data sequence around the source, the times in minutes since 01:00 UT are 
given in the title, location of time-distance sources is marked by orange star in right-most 
panels. Middle row: HMI magnetogram (left) and velocity images with the lines along which 
the stack plots are made. The values of and 1 correspond to horizontal lines (lower and upper 
correspondingly) in the middle row plots; The red arcs represent angle of integration used for 
computing time-distance diagrams in Figure (2] Red stars mark locations of time-distance 
sources. Bottom row: stack plot from velocity data. In all velocity plots, black and white 
correspond to downward and upward motions respectively. The black in magnetogram (middle 
row, left) corresponds to negative polarity, white to positive. 
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acoustic energy estimates in Section[3]we use the same egression computation but 
with 1-mHz frequency bandwidth) and using Green functions built for surface 
monochromatic point source of corresponding frequency using geometrical optics 



approach (Donea, Braun, and Lindsey, 1999 Donea, Lindsey, and Braun, 2000 



Lindsey and Braun, 2000] Matthews, Zharkov, and Zharkova, 2011 



Zharkov, 



Green, Matthews et a/., 2011; [Zharkov, Zharkova, and Matthews, 2011). The 



pupil size is set from 10 to 40 Megameters. As flare acoustic signatures can 
be submerged by ambient noise for the relatively long periods over which the 
egression power maps are integrated, we follow [Donea and Lindsey] (|2005p and 
use egression power "snapshots" to discriminate flare emission from the noise. 
Such a snapshot is simply a sample of the egression power within a time At = 



1 



500 seconds. 



2 mHz :— — II 

Time-distance diagrams ([Kosovichev and Zharkova] I1998| IKosovichev1l2006^ 
[20071 IZharkova and Zharkovi [20071 [K osoviche v] [20111 Zharkov. Zharkova, and 
Matthews H^^^ are computed bv selecting a source location, rewriting the 
observed surface velocity signal [v] in polar coordinates relative to the source 
[v{r^0^t)] and then using the azimuthal transformation 

Vm{r,t)= / v{r,0,t)e-'^'dO, (3) 

J a 

where integration is normally performed over the whole circle (i.e. a = 0, /3 = 
27r). In the case when m = 0, the integration limits can be chosen to represent an 
arc. Acoustic wave-packets of a sufficiently strong amplitude relative to the noise 
are seen as a time-distance ridge that largely follows a theoretical time-distance 
curve. 

The time-distance method is an observationally direct technique. If one is 
to think of circular ripples propagating from a source, at any moment of time, 
the procedure basically contracts ripples at a particular distance into a point 
on a time-distance diagram. Since only a handful of known sunquakes were 
accompanied by visible ripples, this obviously increases the chances of detection. 
Nonetheless, the detection of the time-distance ridge can be affected by many 
factors such as noise in the data, filtering and image-processing techniques used. 
The acoustic egression method, on the other hand, provides a more quantitative 
measurement of the seismic source, but relies on a theoretical model of acoustic 
wave propagation from the source through solar interior to the surface. In a way, 
this method contracts the time-distance ridge into a single point measurement. 
Indeed, if we imagine a point source momentarily generating acoustic waves in 
an ideally quiet Sun, e.g. without convective motions, then these waves will be 
observed as accelerating ripples on the surface, contracted to a ridge in a time- 
distance diagram computed at the source, and seen as a bright emission around 
the source time and location in an egression power map. The holographic method 
is more sensitive, but vulnerable to noise, variations from the assumed model and 
is susceptible to an increased possibility of false detections, which is normally 
handled by analysing the statistical significance of the detected signal ( Donea, 
Lindsey, and B raun, 2000; [Donea and Lindsey, 2005[ Matthews, Zharkov, and 
Zharkova, 201 it 
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Figure 4. Frequency response for Source 1 (left) and Source 2 (right). Egression power rms as 
function of time (along x-axis) and frequency (along y-axis) are plotted in units of quiet-Sun 
egression power (top row) and as log^o of acoustic energy flux [ergcm~^ s"-*^] (bottom row) 
for corresponding frequency band. 

3. Results 

Figure [TJa) shows the egression power snapshot in the 6-mHz frequency band at 
around 01:51 UT. The data is saturated at three times the quiet-Sun egression 
power for better contrast. The image shows the locations of the two seismic 
sources as bright regions. The larger source (on the left, East) corresponds to the 
one reported by IKosovichevI (j2Q 1 1 p . We will call it Source 1. The smaller source 
on the right (West, Source 2) was detected by |Zharkov, Green, Matthews et al 
(|2Qlip . Figures [TJb) - (e) provide context for seismic observations by overplotting 
egression-power contours at 2.5 and 3 times the quiet- Sun egression units on HMI 
intensity, magnetogram and velocity data. 

In Figure [2] we present time-distance diagrams computed from unfiltered HMI 
velocity difference data for both sources. Locations of the time-distance sources 
are marked by red stars in Figure [TJb) - (e). Theoretical time-distance lines fitted 
to the data provide an estimate for the quake start times, suggesting that Source 
2 is initiated at around 01:49:30 UT, about a minute earlier than Source 1. 

Time-distance diagrams for both sources are computed in the direction of the 
surface velocity transient movement with integration in ([3]) performed over 90° 
arc. These are shown in Figure [3] where we present the results of our analysis 
of the the horizontal motions of velocity transients at the quake locations. The 
upper plots are a time sequence of velocity data from small regions encompassing 
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each seismic source, showing strong transient downflows and horizontal move- 
ment at the time-distance sources, which are marked by a star in the rightmost 
panel in each case. The middle row plots show red dashed arcs representing the 
^-integration range used for computing time-distance diagrams in Figure [2] (see 
Equation (|3])). These are overplotted on HMI line-of-sight magnetogram and 
magnetogram-difference images. There we also define stack lines along which we 
present the velocity variation in the bottom row. The time is along the x-axis, 
and horizontal distance is along the y-ax.is. In each case, the point marked as 
zero, corresponds to the lower horizontal line (at ^ = 0) in the bottom row, with 
1 corresponding to the upper horizontal line. 

For Source 1 (left column in Figure [3]) we can see a strong downfiow of 
about 500 — 700 m s~^ starting around 01:50 UT followed by the upfiow, trav- 
elling around 2.5 Mm along the line in around two -three minutes. This gives a 
speed estimate between 13 — 20 km s~^, which is around the value reported by 
IKosovichevI (j201ip . Source 2 (right column in the Figure [3]) shows a persistent 
strong downfiow of around 1 km s~"^ located near the zero point from around 
01:48 UT to 01:52 UT as well downward transient movement along the line. The 
horizontal movement starts about 01:48:30 UT and moves with the estimated 
speed 14 - 22 km s"^ 

To measure the acoustic energy for the seismic sources we have recomputed 
the acoustic egression as described above but applying 1-mHz bandwidth fil- 
tering, instead of 2 mHz, so that the measurements at each frequency band 
do not overlap. We used a 6-mHz egression-power snapshot to define kernels 
for spatial integration by thresholding at the quiet- Sun intensity and selecting 
congruent areas corresponding to sources one and two. Then, at each frequency, 
the egression power scaled by the sound-speed and density values taken from 
Model C dChristensen-Dalsgaard, Dappen, Ajukov et a/., 1996D at around 250 
km above the photosphere was integrated over such kernels and a 30-minute 
period around the peak of HXR emission. The results are presented in Figure 
[5l From this we estimate that the total amount of acoustic energy released by 
Source 1 is around 1.18 x 10^^ ergs and Source 2 is around 6.08 x 10^'' ergs. 

Figure |4] shows the rms variation of egression power spatially integrated 
around the two seismic sources as a function of frequency and time. In the top 
row of Figure m the egression power is expressed in units of quiet-Sun egression 
power computed for the corresponding frequency, with the log of the integrated 
egression in ergcm"^ presented in the bottom row. The relative strength of 
the eastern source (Source 1, left column) is apparent, with a significant increase 
in relative egression power seen in all frequencies. The western source (Source 2, 
right column) has a clear signature around 4-7 mHz range, with little apparent 
power increase around 3 mHz or frequencies above 9 mHz. The egression method 
does not provide an accurate determination of the quake onset time due to the 
frequency-bandwidth filtering of the data and the induced timing uncertainty of 
500 or 1000 secon ds for 1-mHz and 2-mHz filters respectively ( Donea, Braun, 
and Lindsey, 1999; |Donea and Lindsey, 2005D . Nonetheless, it is interesting to 
note from Figure H] the relative frequency variation of egression peak times. 

Figure [6] shows the photospheric changes in HMI line-of-sight velocity, mag- 
netic field, and continuum intensity at the time-distance sources (indicated in 
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Figure 5. Acoustic energy distribution versus frequency for the two sources. Solid line 
corresponds to Source 1, dashed to Source 2. 

Figure [1]) over two hours (top panel) and 40 minutes (bottom panel) around the 
flare. 



4. Discussion 

When magnetic field is present, as is the case for sunspots, the acoustic waves 
become magneto- acoustic. Numerical and theoretical modelling of acoustic- wave 
packet pro pagation through a magnetised plasma (jCally, 2QQQ[ Schunker and 
Call y, 2QQ6;|Shelyag, Zharkov, Fedun et al, 2QQ9| Moradi, Baldner, Birch et al, 
2010; [Felipe, Khomenko, and Collados, 2Q1Q[ for example) shows that acoustic 
waves when passing through magnetic-flux tube transform into different types 
of magneto- acoustic waves, which we then expect to flnd in sunspots. Mode con- 
version between various magneto-acoustic waves takes place around the region 
where the ambient sound speed [c] equals the Alfven velocity [a] which is situated 
near the surface in sunspots. For instance, the slow-mode waves (classifled as 
being slow below the a = c transition region) have properties fundamentally 
different from those expected from non-magnetic acoustic waves. These waves 
propagate along magnetic-fleld lines and therefore have very different wavefronts 
and speeds. However, as shown by |Shelyag, Zharkov, Fedun et a/.|(|2009p . the fast- 
mode magneto- acoustic waves are similar to acoustic waves. While the presence 
of a magnetic fleld and associated changes in ambient conditions clearly affect 
their paths and speed, these changes are relatively small. When travelling outside 
the magnetic tube, these waves become acoustic. 
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While the physical mechanism behind the excitation of the waves remains un- 
determined, we can hypothesise that the observed flare-generated acoustic wave- 
fronts are likely connected to such magneto-acoustic waves. Indeed, as quakes 
happen in magnetised plasma, the generated waves will be magneto- acoustic. 
However, the acoustic nature of the ripples and ridges in the time-distance dia- 
grams is very clear: detected ridges closely follow the theoretical travel time (see 
IKosovichev I (|2QQ6l |2QQ71 [2QI1|; .Martmez-Qliveros, Moradi, and Donea] (EOUHl), 
e.g.). Nonetheless, the subsurface conditions will affect the time and speed of the 
propagation of such waves, introducing an anisotropy of the wave front such as 
that reported bv IKosovichevI (j2Q 1 1 p for Source 1. As the egression measurements 
are based on a quiet-Sun model of acoustic-wave propagation, they are certainly 
affected by such an anisotropy. On the other hand, the travel time perturbations 
of acoustic waves passing through sunspots measured by time-distance helio- 
seismology ([Duvall, Jefferies, Harvey et al.^ 1993[ Duv all, Kosovichev, Scherrer 
et a/., 1997; [Kosovichev, Duvall, and Scherrer, 2000; Zharkov, Nicholas, and 
Thompson, 2007) indicate that such changes are expected to be small and of the 
order of one minute. 

The stacked nature of egression kernels presented in Figureflja), is very similar 
to t he measurements obtained for other quakes ( Donea, Braun, and Lindsey, 
1999; [Donea and Lindsey, 2005[ [Donea, 2011D . In these works the authors suggest 
that this could be the result of interference caused by the rapid motion of the 
source, roughly in the direction along which the kernels are stacked. For example, 
Donea and Lindsey] (|2005p found that for 28 and 29 October 2003 quakes the 



motion of the HXR sources was indeed aligned accordingly with egression power 
stacks. This was also confirmed for the 23 July 2002 flare bv .Kosovichev ( 20Q7|) 
using HXR and Doppler data. 

For this flare, 15 February 2011. IKosovichevI (|201ip used time-distance dia- 
gram analysis to detect a supersonic (adiabatic sound speed at the photospheric 
level is around 7 — 8 km s~^ in the quiet Sun) movement of the location of Source 
1 of around 15 — 17 km s~^. In fact, apparent horizontal motions of downward 
transients are observed in the HMI velocity-data sequence at the location of 
both seismic sources before and around the time of the quakes (see Figure [3] for 
instance). In order to investigate this, we have produced stack plots (bottom 
row of Figure [3]), from which the horizontal speed in the direction along the 
integration line is estimated (Section[3]). While the sound speed in the penumbral 
magnetised plasma is likely to be different from the quiet Sun, given the estimates 
(13 - 20 and 14 - 22 km s"^ for sources 1 and 2 respectively) it is reasonable to 
conclude that in both cases the motion appears to be supersonic. We note that 
supersonic horizontal movements in the photosphere have also been reported in 
association with several other flare induced sunquakes (see for example Donea 
and Lindsev n2Qb5[ IKosovichevI I2007p . 

Considerable anisotropy in the acoustic amplitude of the ripples from the 



vantage of the sources has been observed for most quakes (Kosovichev, 2006 



Moradi, Donea, Lindsey et al.^ 2007|[Donea, 201ip . In fact . IDoneal (|20 1 1 p suggests 



that the maximum amplitude of the ripples emanating from a moving source is 
generally along the axis of the source, displaced from the source location in 
the direction of the motion. This is confirmed for both of our sources by our 
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Figure 6. Variations of velocity, intensity, and magnetic field at time-distance sources 1 (left) 
and 2 (right); Top plots show two-hour time series around the time of the flare, the bottom 
plots are the same but on a 45 minute scale. The red dotted line in the top plots indicate the 
edges of the bottom plot for reference. The vertical dashed and solid blue lines correspond to 
01:49:57 and 01:55:57 UT respectively. Please see text for more details. 



directional time-distance diagrams (see Figures [2] and [3l where after considering 
various arcs, the results where the time-distance ridge appears the strongest are 
presented). 

As mentioned above, theoretical time-distance curves fitted to the ridge data 
in the time-distance diagrams suggest a marginally earlier start for the west- 
ern quake. As this difference is within a margin of error, we consider velocity 
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plots from both time-distance sources shown in Figure [6l These indicate strong 
downflows around the times of the quake. It is also clear that the Source 2 
downflow happens about 45 seconds earlier than Source 1. This confirms the 
timing difference between the two sources. 

The plots for the magnetic-field and intensity variation over two hours for 
Source 2 in Figure [6] show a long-term continuous increase in magnetic-field flux 
density coupled with the intensity decrease, which indicates ongoing emergence 
of the flux at the location. Around the time of the flare, however, we see a gradual 
and apparently permanent magnetic-field increase of around 130 — 160 Gauss, 
taking place from around 01:50 to 01:56 UT, coinciding with the peak of the 
flare X-ray emission. This is preceded by a 7% increase in intensity taking place 
approximately from 01:47 to 01:50 UT. After its peak, the intensity decreases 
slowly to a pre-flare level, suggesting that a physical process not related to sudden 
heating, but perhaps associated with the magnetic-field change, is taking place. 
There is also present a transient magnetic-field variation of about 50 Gauss at 
around 01:51-01:52 that is likely to be associated with the abrupt changes in 
the observed line profile due to precipitation of energetic particles. An abrupt 
increase in the line-of-sight magnetic field could indicate that flux-rope field 
lines have become more vertical, which would be in agreement with the eruption 
model suggested by [Zharkov, Green, Matthews et aL\ (|201ip . Source 1 plots, 
on the other hand show a transient increase in the magnetic-field strength of 
about 50-60 Gauss around 01:50 followed by a rapid decrease with the field 
settling back to pre-flare levels about half an hour later. There is also a very 
small increase in continuum intensity increase of around 2.7%. 

Similar plots of photospheric variations where the data were averaged over 
the egression so urces have also been presented in Figure 4 of Zharkov, Green, 
Matthews et al. (|201ip . We note the difference between the two figures in the 
magnetic-field and intensity responses. We find that these parameters vary sig- 
nificantly from pixel to pixel in egression kernels, suggesting that it is very 
important to pinpoint the exact location of the quakes. These can be determined 
more precisely using time-distance diagrams, as pixels with respect to which 
the diagrams showing a clear ridge are computed form a connected set for each 
source. However, we see different magnetic-field behaviour at the seismic location 
from pixel to pixel: for example, a few pixels show signs of a permanent magnetic 
change, many show an oscillatory kind of behaviour, and even more show a 
transient-like response. Indeed, the results of our preliminary analysis suggest 
that abrupt and permanent magnetic changes are more prevalent in the ribbons, 
although the picture is still quite complex, perhaps due to the continuing flux 
emergence in the re gion. Further analysis making use of known methods ( Sudol 
and Harvey, 2005; Z harkova, Zharkov, Ipson et al.^ 2005D and now available HMI 
vector magnetic data is needed to understand the variability of magnetic-field 
changes in seismic sources and this flare in general. 

The acoustic-energy estimates obtained here are somewhat higher than ob- 
taine d for the same flare by Alvar ado- Gomez, Buitrago-Casas, Martmez-Oliveros 
et fl/.' |2M2| . This is most likelv due to several factors such as use of Model C as 
opposed to VAL, different pupil dimensions and different integration kernels (we 
use larger areas around 37 and 33 Mm^ for Sources 1 and 2 correspondingly, as 
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opposed to 12 Mm^ ( Alvar ado- Gomez, Buitrago-Casas, Martinez- Oliveros et al. 
(|2012p , C. Lindsey, private communication, 2012). The estimate for Source 1 puts 
it amongst the most powerful sunquakes associated with X-class flares (Moradi, 
Donea, Lindsey et al, 2007; [Besliu-Ionescu, Donea, Cally et a/., 2QQ5D . We note, 
however, that [Besliu-Ionescu, Donea, Cally et a/.|(j2006|) reports an even stronger 
seismic event for the M7.8 flare on 2 December 2005. 



5. Conclusions 



We have presented here a time-distance diagram in addition to the one found 
by IKosovichevI (j201ip showing a clear ridge for the second seismic source associ- 
ated with the 15 February 2011 X-class flare. Using time-distance analysis and 
HMI line-of-sight velocity observations we deduce that the quakes are excited 
at around 01:50 UT, with the eastern source onset preceding the western one by 
about 45 — 60 seconds. We have also detected apparent horizontal motions of 
the downward velocity transients at the time and location of both quakes. The 
speed of such motions is larger than the ambient sound speed. The direction 
of such motions is aligned with the stacked egression kernels and amplitude 
anisotropy of the generated wavefront, indicating that a moving source is the 
likely scenario for both quakes. We estimate the acoustic energy released by 
both quakes to be around 1.18 x 10^^ ergs for Source 1 and 6.08 x 10^^ ergs 
for Source 2. For Source 1 this is about an order of magnitude higher than the 
Lorent z-force energy estimate for a generic flare provided by Hudson, Fisher, and 
Welsch ^2008p . This is i n line with fl ndings by Alvarado- Gomez , Buitr ago- Casas , 
Martmez-Oliveros et al. (|2012p , where more accurate evaluation of Lorentz-force 
energy has been produced. However, given a number of simpliflcations used in 
obtaining the Lorentz-force estimate, such as the use of line-of-sight magnetic 
field only, the assumption of a single area where changes occur and relatively 
low magnetic-field strength, in our view it would be premature to discard the 
Lorentz force as a possible mechanism for quake excitation with further analysis 
based on ( [Fisher, Welsch, and Abbett, 2012] [Fisher, Bercik, Welsch et a/., 2012^ 



making use of fill vector magnetic-field data necessary. 

Further analysis is required in order to understand the physical nature of 
both detected quakes and their link with the fiux-rope eruption that was asso- 
ciated with the X-class flare. In particular, the HMI vector magnetogram data 
should shed the light on full magnetic changes, and numerical extrapolations 
of the three-dimensional magnetic fleld from photosphere through atmosphere 
and corona will give us a clearer picture of the magnetic-fleld restructuring and 
energy release associated with this event. 
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